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SUlOfARY 


Approximate  analytical  relations  are  developed 
for  several  possible  "conical"  flow  fields  resulting  from 
the  steady  supersonic  flow  of  a  uniform,  non-vlscous,  non¬ 
heat  conducting,  combustible-gas  mixture  past  a  semi- 
Inflnlte  cone  (at  zero  angle  of  attack)  for  Instantaneous 
chemical  reaction.  The  flow  fields  considered  are  of  the 
shock-deflagration,  detonation,  and  detonation-shock  types 
with  attached  adiabatic  and  diabatic  discontinuities. 
Particular  emphasis  Is  placed  upon  the  Chapman- Jouguet 
detonation. 


Parametric  curves  are  presented  for  shock- 
deflagration  flows  about  cones  and  wedges,  and  Chapman- 
Jouguet  detonatlve  flows  about  cones. 

Generalized  analyses  of  the  flow  characteristics 
of  an  oblique  exothermic  discontinuity  are  given  in  the 
appendices . 
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SYMBOLS 

local  speed  of  sound 
normalized  local  speed  of  sound,  a/c 
maximum  speed  obtainable  by  adiabatically 
expanding  into  a  vacuum  (constant  for  any 
given  region  0,1,2,  etc.)  c  ■  \l  2jrRr/()'-l) 
constants  defined  following  Eq.  (7) 

functions  defined  respectively  by  Eqs.  (B4) 
and  (B6) 

specific  heat  at  constant  pressure 
specific  heat  at  constant  volume 
cot 

s 

cot  <Pa 

local  Mach  number,  V/a 

component  of  local  Mach  number  normal  to  a 
conical  ray,  v/a  (see  Fig.  2) 
static  pressure 
total  pressure 

mechanical  equivalent  of  heat  added  per  unit 
mass  of  fluid  (ft-lb/slug,  or  other  consistent 
units) 

heat-addition  parameter  (Damkohler's  second 
parameter) 
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SYMBOLS  (coDt'd) 

gas  constant 
static  temperature 

velocity  components  parallel  and  perpendicular 
respectively  to  a  conical  ray  (see  Fig.  2) 
normalized  velocity  components,  u/c,  v/c 

/  2  2 

local  resultant  velocity,  V  u  +  v 
free-stream  velocity 
(r  -  l)/2 

ratio  of  specific  heats,  c  /c 

p  V 

local  flow  inclination  relative  to  free- 
stream  direction  (see  Fig.  2) 
mass  density 
total  temperature 

spherical  or  polar  angle  (see  Fig.  2) 
the  maximum  body  semi-apex  angle  for  which 
a  conical  flow  region  can  exist  for  a  given 
free-stream  Mach  number  and  heat  release 
body  semi-apex  angle  corresponding  to  a 
surface  Mach  number  of  unity 
functions  of  y  defined  in  Appendix  C 
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0,1, 2, etc. 


B 

D 

d 

F 

J 


denotes  a  constant-value  quantity  in  the 
particular  region  (see  Fig.  1)  indicated 
by  the  subscript 

denotes  conditions  at  the  body  surface 
denotes  conditions  at  a  detonation  wave 
denotes  conditions  at  any  general  exothermic 
discontinuity 

denotes  conditions  at  a  deflagration  wave 
denotes  conditions  at,  or  pertaining  to,  a 
Chapman- Jouguet  detonation 

denotes  conditions  at  a  singular  discontinuity 
denotes  conditions  at  an  adiabatic  shock  wave 
denotes  juncture  point  of  inner  and  outer 
solutions  for  the  unit  field 
see  statement  following  Eq.  (Bib) 

see  statement  following  Eq.  (B6) 

Superscripts 


indicates 
(see  Fig. 
indicates 
(see  Fig. 


a  variable  quantity  in  region  0 
1) 

a  variable  quantity  in  region  1 
1) 
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Superscripts  (cont'd) 
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=  Indicates  a  variable  quantity  in  region  2 

(see  Fig.  1) 

(J)  denotes  conditions  corresponding  to  a 

Chapman- Jouguet  detonation 
(o)  see  statement  proceeding  Eq .  (B9) 


5 


Ttw  JohM  Hafklrw  Unhwtliy 
>ru»  FHVWM  lAM««T«aV 
Ilivtr  Sprint,  ttarylwt 


INTRODUCTION 

Although  supersonic  combustion  In  the  form  of 
a  non-statlonary  Chapman- Jouguet  detonation  propagating 
through  a  tube  filled  with  a  combustible  gas  has  been 
studied  by  combustion  scientists  for  about  eighty  years 
(Refs.  1  and  2),  It  Is  only  In  recent  times  that  experi¬ 
mental  Investigations  have  been  conducted  on  supersonic 
combustion  processes  which  are  stationary  with  respect 
to  an  aerodynamic  body,  a  fuel  Injection  nozzle,  or  an 
exhaust  nozzle  (see,  for  example.  Refs.  3  through  8). 

The  motivation  for  the  Increased  research  activity  In 
supersonic  combustion  stems  from  recent  analyses  of  hyper¬ 
sonic  ramjets  utilizing  supersonic  combustions  (Refs.  9, 

10,  and  11),  proposals  for  detonatlve  propulsion  engines 
(Refs.  12  and  13),  and  studies  of  unconventional  supersonic 
aircraft  utilizing  external  burning  for  lift  and  for 
propulsion  (Refs.  14  and  15),  Hence,  there  are  now 
practical  as  well  as  fundamental  reasons  for  obtaining  a 
more  complete  understanding  of  supersonic  combustion 
processes , 

It  is  observed,  with  regard  to  the  detonatlve 
supersonic  combustion  experiments  of  Refs.  4,  5,  and  6  that 
blunt  bodies  are  more  conducive  to  Initiating  shock-induced 
combustion  than  pointed  ones  because  of  the  relatively  higher 
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Static  temperature  behind  the  normal  shock.  On  the  other 
hand,  It  seems  only  logical  that  greater  dividends  are 
derived  from  experimental  Investigations  performed  on  con¬ 
figurations  amenable  to  analytical  calculation.  A  wedge 
Is  much  more  satisfactory  In  this  respect  than  a  blunt  body. 
A  cone,  however.  Is  not  only  amenable  to  analytical  calcula¬ 
tion  (at  least  for  the  Idealized  case  of  Instantaneous 
chemical  reaction)  but  also  represents  a  configuration  which 
Is  adaptable  either  to  the  technique  of  firing  In  a  gas- 
filled  ballistic  range  or  to  stationary  testing  in  a  com¬ 
bustion  wind  tunnel.  This  overlapping  of  available  testing 
techniques  can  hardly  fail  to  be  rewarding.  Also,  a  cone 
is  free  of  the  end  effects  encountered  on  a  finite-span 
wedge.  A  wedge  has  the  advantage,  however,  of  producing  a 
stronger  shock  than  a  cone  for  a  given  apex  angle  and  free- 
stream  Mach  number.  Experimentally  this  can  be  compensated 
for  in  cone  tests  by  increasing  the  temperature  of  the 
ambient  gas  in  the  tunnel  or  range. 

This  paper  presents  approximate  analytical  rela¬ 
tions  for  several  possible  "conical"  flow  fields  resulting 
from  the  steady  supersonic  flow  of  a  uniform,  non-viscous, 
non-heat  conducting,  combustible-gas  mixture  past  a  semi- 
infinite  cone  (at  zero  angle  of  attack)  for  Instantaneous 
chemical  reaction.  Aside  from  its  fundamental  aspects,  the 
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primary  purpose  of  the  paper  is  to  present  material  which 
it  is  hoped  will  be  useful  in  the  design  and  interpreta¬ 
tion  of  cone-flow  supersonic  combustion  experiments  aimed 
at  obtaining  a  better  understanding  of  the  basic  phenomena 
involved  in  supersonic  combustion. 

It  appears  that  the  only  other  theoretical  work 
dealing  with  non-linear  diabatlc  cone-flow  fields  is  that 
of  Kvashnina  and  Chernyi  (Ref.  16).  Kvashnina  and  Chernyi, 
without  recourse  to  numerical  calculation,  deduce  from  an 
examination  of  the  hodograph  differential  equation,  the 
detonation  hodograph,  and  the  cone  boundary  condition,  the 
general  character  of  the  flow  field  for  the  flow  of  a  super¬ 
sonic  gas  stream  past  a  solid  cone  with  an  attached  detona¬ 
tion  wave  and  varying  amounts  of  heat  release. 

Although  supersonic  combustion  flow  about  a  cone 
is  the  main  topic  of  this  paper,  some  numerical  results  for 
wedge  combustion  flow  fields  are  included  as  a  matter  of 
interest . 


GENERAL  DESCRIPTION  OF  SOME  POSSIBLE 
CONICAL  FLOW  REGIMES  FOR  INSTANTANEOUS 
CHEMICAL  REACTION 

The  assumption  of  Instantaneous  chemical  reaction 
allows  a  combustion  zone  (detonation  or  deflagration)  to  be 
treated  as  an  exothermic  surface  of  Inflnltesmal  thickness 
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across  which  the  flow  properties  and  chemical  compositions 
are  discontinuous.  An  Illustration  showing  some  possible 
theoretical  supersonic  conical  flow  fields  with  an  exothermic 
discontinuity  appears  In  Fig.  1.  The  detailed  notation, 
applied  to  a  particular  case  of  Fig.  1,  Is  Illustrated  In 
Fig.  2.  The  adiabatic  shocks  appearing  In  Fig.  1  are 
assumed  not  to  produce  a  change  In  the  specific  heats  or 
In  chemical  composition.  The  different  chemical  compositions 
In  the  regions  upstream  and  downstream  of  an  exothermic  dis¬ 
continuity  are  frozen  throughout  their  respective  fields  and 
the  gas  Is  assumed  to  be  calorlcally  and  thermally  perfect 
throughout  a  given  region  although  It  experiences  an  Instan¬ 
taneous  change  In  the  specific  heats  In  crossing  an  exothermic 
discontinuity.  The  bounded  flow  fields  between  discontinuities 
(shocks,  deflagrations,  detonations,  and  bodies)  are  conical 
In  the  sense  that  the  flow  properties  are  constant  along 
straight  line  elements  emanating  from  the  apex  of  the  body. 

This  is  representative  of  the  flow  field  for  a  chemically 
reacting  mixture  with  finite  reaction  time  in  the  region  far 
away  from  the  apex.  In  the  vicinity  of  the  apex,  for 
distances  of  the  order  of  the  chemical  relaxation  distance, 
the  flow  will  not  be  conical.  The  reduction  of  the  flow 
problem  to  one  Involving  conical  fields  represents  a  con¬ 
siderable  simplification  in  the  analytical  task. 
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For  an  ambient  stream  of  a  given  chemical  com¬ 
position  the  consequences  of  the  chemical  reactions  are 
Implicitly  Included  in  the  aerothermodynamlc  parameters 
of  the  oblique  exothermic  discontinuity  relations  given 
In  Appendix  A.  Since  it  Is  beyond  the  scope  of  this 
paper  to  treat  the  thermochemical  aspects  of  the  flow  In 
detail,  those  quantities  which  would  be  determinate  If  the 
governing  thermochemical  relations  were  Included  are 
assumed  to  be  known  "a  priori".  These  Include,  the  heat 
release,  the  specific  heats  downstream,  and  the  normal- 
component  Mach  number  upstream  of  the  exothermic  dis¬ 
continuity  . 

The  oblique  exothermic  discontinuities  shown  in 
Fig.  1  are  classified  by  applying  Jouguet's  rule  (Ref.  17) 
for  one-dimensional  reactive  flow  to  the  normal  components 
of  the  upstream  and  downstream  Mach  numbers  at  the  dis¬ 
continuity.  Hence,  immediately  ahead  of  ah  oblique  exothermic 
discontinuity  the  normal  Mach-number  component  relative  to 
the  discontinuity  is  subsonic  for  a  deflagration  and  super¬ 
sonic  for  a  detonation.  The  downstream  normal  Mach-number 
component  is  subsonic  for  a  weak  deflagration  or  a  strong 
detonation,  and  supersonic  for  a  strong  deflagration  or  a 
weak  detonation.  Oblique  Chapman-Jouguet  deflagrations  and 
detonations  are  those  for  which  the  downstream  normal  Mach- 
number  component  is  unity. 
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The  use  of  the  adjectives  "strong”  and  "weak"  in 
the  foregoing  classification  leads  to  an  ambiguity  in  termi¬ 
nology  if  the  traditional  method  of  describing  an  oblique 
adiabatic  shock  due  to  a  cone  or  wedge  as  strong  or  weak  in 
ac^|ii?^nce  with  the  magnitude  of  its  Inclination  angle  is 
carried  over  to  the  dlabatlc  shock  case,  since,  in  a  manner 
analogous  to  adiabatic  flow,  two  detonation-wave  angles  are 
obtained  for  a  cone  or  wedge  at  a  given  free-stream  Mach 
number,  and  a  given  detonative-wave  heat  release.  The 
ambiguity^  Is  avoided  by  designating  the  detonation  wave 
with  the  larger  inclination  as  "strongly  inclined"  and  the 
one  with  the  lesser  inclination  as  "weakly  inclined".  Hence, 
an  oblique  strong  detonation  may  be  a  weakly- inclined  strong 

I 

detonation  or  a  strongly-inc^ned  strong  detonation.  Similar 
descriptions  are  applied  to  weak  detonations  and  to  C-J 
detonat ions . 

It  can  be  shown  from  entropy  considerations 
(Ref.  17)  that  strong  deflagrations  are  impossible.  For  one¬ 
dimensional  detonations,  the  literature  is  controversial 
(Refs.  1  and  2)  regarding  the  theoretical  existence  of  non 
C-J  detonations,  although  the  preponderant  opinion  seems 


^This  situation  does  not  occur  for  oblique  adiabatic  shocks 
with  constant  specific  heats,  since  consideration  is  never 
given  to  classifying  such  shocks  in  terms  of  their  normal- 
flow  properties  because  the  normal-flow  weak  solution  yields 
the  result  that  "nothing  happens"  see  Eq .  (A6) . 


11 


Tlw  Mm  HapklM  UnKmilty 
«MUII  rNVaiM  k*MIAT«aV 
SMvw  i#ring,  Maryland 


to  be  that  Chapman -Jouguet  detonations  are  the  most  probable, 
and  that  strong  detonations  are  more  probable  than  weak  ones. 
For  plane  (wedge)  flow,  Slestrunch  and  associates  (Ref.  18) 
have  theoretically  demonstrated  the  Impossibility  of  weak 
detonations  from  an  examination  of  the  detonation  polar 
without  recourse  to  entropy  considerations.  For  cylindrical 
(cone)  flow  Kvashnina  and  Cherny!  (Ref.  16)  imply  the  non¬ 
existence  of  weak  detonations.  In  any  case,  the  present 
analysis  Is  generally  confined  to  weak  deflagrations,  strong 
detonations,  and  Chapman- Jouguet  detonations  all  having  weak 
Inclination  angles.  The  restriction  to  weak  inclinations 
follows  by  analogy  to  oblique  adiabatic  shocks  for  which  only 
weak  Inclinations  are  observed  in  physical  reality. 

For  plane  flow,  the  shock-deflagration  and  the 
detonation  flow  regimes  shown  in  Fig.  1  have  been  analyzed 
previously  in  Refs.  19  and  20,  and  Refs.  16,  21,  and  22 
respectively.  A  relatively  superficial  examination  of  the 
problem  is  sufficient  to  convince  one's  self  that  there 
also  exist  cylindrical  flows  having  the  same  general 
character,  although  the  field  details  differ  in  that  for 
cylindrical  flows,  the  pressure  increases  in  the  downstream 
direction  along  curved  streamlines  in  bounded  conical-flow 
regions,  whereas  for  plane  flows  the  pressure  is  constant 
on  straight  streamlines  in  bounded  conical-flow  regions. 
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For  the  detonatlve  type  flows  it  is  pertinent  to 

note  that  there  are  two  particular  values  of  the  body  semi- 

apex  angle  cp^  (for  both  cones  and  wedges)  which  are  crucial 

to  the  delineation  of  the  detonatlve-f low  regimes  shown  in 

Fig.  1.  For  given  values  of  the  free-stream  Mach  number, 

specific  heats,  and  heat  release  at  the  detonation  wave, 

there  exists  a  maximum  body  semi-apex  angle  (Pq  above 

max 

which  a  conical  flow  (plane  or  cylindrical)  cannot  exist, 

l.e.,  the  wave  becomes  detached,  and  there  exists  also  a 

semi-apex  angle  qPo  ,  cPn  >  'Po  >  which  the  detonation 
®  J  max  ~  ®  J 

wave  becomes  a  Chapman- Jouguet  detonation.  For  body  semi- 

apex  angles  less  than  but  greater  than  (pg  the  flow 

max  **J 

equations  yield  two  solutions,  as  previously  noted,  one  of 

which  corresponds  to  a  strongly-inclined  detonation  wave  and 

the  other  to  a  weakly-inclined  detonation  wave.  For  body 

semi-apex  angles  less  than  (pQ  ,  the  detonative  wave  angle 

®J 

remains  at  qPj,  and  an  expansion  region  occurs  downstream  of 
the  detonation  wave.  This  expansion  region  has  a  different 
character  for  plane  and  cylindrical  flows.  For  plane  flow, 
the  expansion  region  is  simply  a  Prandtl-Meyer  expansion  of 
sufficient  extent  to  turn  the  flow  parallel  to  the  body 
surface.  For  cylindrical  flow,  Kvashnina  and  Chernyi  in 
Ref.  16  ascertain  that  the  conical  expansion  field  adjacent 


13 


Tlw  Mmm  Uffhmtlly 

«^UM  riivwM  UMMMTtav 
Stivar  iprlnf,  M(rvltn4 


to  the  detonation  wave  Is  bounded  on  the  downstream  side 
by  a  shock  wave  followed  by  a  compression  flow  field  between 
the  shock  and  cone  (Fig.  Id).  As  the  cone  angle  decreases, 
the  width  of  the  expansion  zone  Increases,  while  the  In¬ 
tensity  of  the  shock  wave  at  first  Increases  and  then  begins 
to  decrease  until  It  degenerates  Into  a  characteristic. 

The  flow  direction  on  this  characteristic  Is  parallel  to 
the  axis  of  symmetry.  An  Isentroplc  compression  then 
occurs  between  the  characteristic  and  the  cone  surface. 

When  the  cone  angle  is  decreased  still  further  the  amount 
of  compression  Is  reduced  until  at  zero  cone  angle  It  dis¬ 
appears  entirely  and  the  region  behind  the  conical  expansion 
region  is  uniform  and  parallel  to  the  axis  of  symmetry. 

For  a  cone,  the  detonative  flow  in  Fig.  Ic  bears 
a  close  relationship  to  the  adiabatic  flow  of  Fig.  la  as 
now  will  be  demonstrated. 

If,  In  an  adiabatic  flow  with  a  given  cone  angle, 
normalized  cone-surface  velocity  Ug,  and  specific  heat  ratio, 
the  normalized  governing  differential  equation  (see  Eq.  (1)) 
is  Integrated  forward  of  the  cone,  ignoring  the  upstream 
boundary  condition,  an  upper  bound  to  the  extent  of  the 
conical-flow  field  is  found  to  occur  at  the  singular 
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discontinuity^  surface,  cp  ,  shown  In  Fig,  la  (Ref.  25). 

s 

This  singular  behavior  In  the  differential  equation  occurs 

at  V  -  a,  which  Is  Identical  to  the  Chapman- Jouguet  condition 

for  dlabatlc  flow.  Therefore,  when  the  detonation  In  Fig.  Ic 

Is  a  C-J  detonation  (<Pq  -•  <Pj),  the  flows  In  the  bounded 

regions  (<Pg^  <P  <Pg)  and  (<Pj <P  ^  <Pg)  of  Figs,  la  and  Ic 

respectively  are  identical  for  a  specified  cone  angle, 

normalized  cone-surface  velocity,  and  specific  heat  ratio. 

It  Is  well  known  (Ref.  23)  that  for  adiabatic  cone  flow 

there  exists  no  uniform  upstream  flow  that  can  be  matched 

to  the  conical  flow  region  bounded  by  the  solid  cone  and 

the  singular  discontinuity  surface.  However,  for  a  dlabatlc 

flow  with  known  values  of  the  specific  heats  on  both  sides 

of  the  exothermic  discontinuity,  it  is  possible,  for  cone 

angles  less  than  cp-  ,  to  find  a  uniform  upstream  flow 

“max 

having  a  Mach  number  and  total  temperature  compatible  with 
the  existence  of  a  Chapman- Jouguet  exothermic  discontinuity 
at  the  singular  surface.  Furthermore  for  any  assumed  wave 
angle  <p„,  less  than  cp ,  but  greater  than  cp  ,  It  Is  possible 
to  match  that  portion  of  the  foregoing  conical  region 


Also  called  the  "limit  line"  (Ref.  23)  or  "limit  cone" 
(Ref.  24).  Ref.  25  calls  this  discontinuity  an  "outer" 
discontinuity,  since  there  exists  also  an  "inner"  dis¬ 
continuity  in  the  interior  of  the  cone.  This  distinction 
will  not  be  made  herein  since  the  "inner" discontinuity 
is  not  pertinent  to  the  present  problem. 
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bounded  by  and  cp^  to  a  uniform  free  stream  having  a 
Mach  number  and  total  temperature  compatible  with  an  oblique 
detonation  of  Inclination  cp^. 

In  view  of  the  significant  role  played  by  the 
conical  field  bounded  by  a  singular  discontinuity  and  a 
solid  cone  It  Is  convenient  to  refer  to  such  a  field  as  a 
unit  conical  field,  or  unit  field  (since  v/a  varies  from 
0  to  1.0).  A  unit  field  is  uniquely  defined  by  specifying 
<Pg,  Ug,  and  or  cpg,  u*,  and  y- 

One  of  the  purposes  of  this  paper  is  to  derive 
approxlhiate  relations  for  the  unit  conical  field  from  which 
the  detonative  flows  of  Fig.  Ic  can  be  obtained.  This  has 
been  done  by  developing  a  solution  valid  in  the  region  of 
the  singular  discontinuity  and  one  valid  In  the  region  of 
the  cone  surface  and  joining  the  two  together  at  an  inter¬ 
mediate  location.  The  solution  valid  near  the  singular 
discontinuity  provides  also  the  initial  field  Just  downstream 
of  the  C-J  detonation  in  the  detonation-shock  flow  shown  in 
Fig.  Id. 

The  approach  taken  In  obtaining  the  dlabatlc 
conical  flow  fields  considered  herein  is  as  follows. 

For  the  shock-deflagration  case,  the  free-stream 
Mach  number  and  the  adiabatic  shock  angle  are  specified. 

The  corresponding  adiabatic  flow  downstream  of  the  shock  is 
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obtained  using  one  of  the  several  approximate  methods  avail¬ 
able  in  the  literature.  The  inclination  of  the  deflagration 
is  then  taken  such  that  the  desired  upstream  normal  Mach 
number  component  is  obtained  ahead  of  the  deflagration. 

For  a  given  heat  release  the  flow  field  downstream  of  the 
deflagration  and  the  corresponding  wedge  or  cone  angle  is 
then  determined  by  the  approximate  method  given  herein  in 
the  next  section. 

For  the  detonation  case,  the  free-stream  Ifach 
number,  all  the  specific  heats,  and  the  C-J  detonation  wave 
angle  are  specified.  The  corresponding  cone  angle  is  then 
determined  by  the  approximate  method  given  herein.  That 
portion  of  the  unit  field  bounded  by  the  C-J  wave  angle  and 
the  adiabatic  shock  angle  then  may  be  used  to  determine  the 
upstream  Mach  number  and  total  temperature  corresponding  to 
various  strong  detonation  wave  angles  for  fixed  values  of 
the  cone  angle,  normalized  cone-surface  velocity,  and  unit- 
field  specific  heat  ratio. 

To  the  best  of  the  writer's  knowledge  the  case  of 
a  strong  detonation  on  a  wedge  (Fig.  Ic)  is  the  only  one  of 
the  reactive  cases  shown  in  Fig.  1  to  have  been  observed 
experimentally  (Ref.  6).  An  isolated  attempt  in  Ref.  4  to 
produce  a  detonatlve  reaction  on  a  cone  was  unsuccessful. 
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It  Is  not  known  whether  the  deflagration  case 
shown  In  Fig.  lb  can  be  produced  In  the  laboratory.  Aside 
from  the  uncertainties  of  Ignition  ^ and  burning,  the  fluid 
dynamic  aspects  of  the  flow  are  somewhat  restrictive  with 
regard  to  laboratory  experiments.  As  will  be  seen  later, 
for  plausible  flame  speeds  and  heat  releases,  the  conical 

2 

flow  downstream  of  the  deflagration  Is  partially  subsonic 
for  a  fairly  wide  range  of  free-stream  Mach  numbers  and  cone 
angles.  Such  a  flow,  although  theoretically  acceptable  for 
a  semi-infinite  body,  cannot  be  attained  In  the  laboratory 
because  of  practical  considerations  limiting  the  test  model 
to  a  finite  length.  The  existence  of  a  subsonic  region  on 
a  finite-length  conical  model  results  In  the  destruction  of 
the  conical  flow  as  a  consequence  of  the  upstream  propaga¬ 
tion  of  downstream  disturbances. 

Assuming  that  a  shock-deflagration  flow  (completely 
supersonic  conical  field)  could  be  experimentally  established, 
it  presumably  could  be  used  to  study  (in  a  wind  tunnel,  say) 
the  transition  from  deflagration  to  detonation  by  systematic 

^An  ignition  system  is  necessary. 

2 

Partially  subsonic  conical-flow  regions  exist  also  for  the 
detonation  and  detonation-shock  cases.  In  the  detonation 
case  (which  is  the  one  of  principal  Interest  here),  however, 
the  part  of  the  flow  spectrum  for  which  this  occurs  is 
extremely  small  (see  Fig.  13)  for  weakly-inclined  detonations. 
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variations  In  the  free-stream  composition  and  aerothermo- 
dynamlc  conditions. 


ANALYSIS 

GOVERNING  DIFFERENTIAL  EQUATION 

The  governing  differential  equation  for  non-linear 
conical  flow,  normalized  with  respect  to  the  maximum  speed, 
c,  is  (Ref.  25), 


1(2  ♦  21  1  ♦  ♦  \  * 

V  -  a  )  u  /d?*^  }  -  a  2  u  -  V  cot  (p  -  u  V  (1) 


where 


V  >■  -  du  /d5p 


(2) 


and 


a  -rll-u  -V 


(3) 


THE  UNIT  CONICAL  FLOW  FIELD 


The  Outer  Solution 


If  it  is  assumed  that  the  right  hand  side  of  Eq.  (1) 

is  not  zero  at  cp  «  (p  ,  an  expansion  for  u  which  will  account 

s 

for  the  singular  behavior^  and  the  boundary  condition  at  <p  is 

s 


*  *  *  *2  „ 

u  -  u„  +  a„  Z  +  C. 
s  s 


3  +  C^  Z^  +  Cg  Z®  +  Cg  Z®  + - 


(4) 


3  ^ 

Terms  to  be  at  least  order  of  Z  in  u  must  be  retained  to 
include  the  singular  behavior. 
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from  which 

-  S  +  I  ^3  2  +  2  z2  +  I  Cg  z3  +  3  Cq  — 


where 


Z  a  +  (cPg 


and, from  Eq.  (3), 


(5) 


(6) 


with 


1 

+  (<P„  -  <P)^  <<  1.0 
s 


(7) 


Substituting  Eqs.  (4)  and  (5)  in  Eq.  (1),  following 
some  long  and  tedious  algebraic  manipulations,  the  relations 
for  the  coefficients  in  Eq.  (4)  are  found  to  be 


C  -  i 
^3 


<"s  -  “  “s> 


tlj  u*  +  Q, 


k 


a 


♦ 

s 


C 


5 


/  *2  *  ♦ 
[^4  %  +  ^5  ^ 


+  Og 


k‘ 
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^6 


+  ng 


+ 


f’s  "s 


2  2 
+  "lO  “  “s  ‘a 


+  Ojj  k 


Relations  for  the  omega  parameters  ^2*  which 

are  functions  of  y,  are  given  In  Appendix  C.  Numerical 
values  are  listed  in  Table  I  for  y  *  1.2,  1.3,  1.4,  and 
1.405^. 

For  a  compression  flow  the  minus  sign  Is  selected 
In  the  expression  for  whereas  for  an  expanding  flow  the 
plus  sign  Is  chosen. 

The  Inner  Solution 


Since  V  Is  zero  at  the  cone  surface,  a  reasonable 
approximation  in  the  vicinity  of  the  cone  is 

(v*/a*)  1.0  (8) 

Utilization  of  the  approximation  of  Eq.  (8)  in  Eq.  (1)  results 
in  the  simplification  of  the  governing  differential  equation 
to 

\i* y dcp^  +  cot  CP  (du'Vd'P)  +  2  u*  -  0  (9) 


1 


Value  used  in  Ref. 


25. 
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When  Eq.  (9)  Is  expressed  In  terms  of  cos  9  as  an  Independent 
variable,  It  may  be  recognized  as  Legendre's  differential 
equation  of  degree  one. 

The  general  solution  to  Eq.  (9)  is 

u*  -  A  cos  <p  +  (B/2)  cos  q)  log  -  B  (10) 

from  which 

V*  -  A  sin  (p  +  (B/2)  sin  cp  log  ~  go|'  +  ®  co't  'P  (H) 

and  the  constants  A  and  B  are  determined  at  the  Juncture  of 
the  Inner  and  outer  solutions. 

Equation  (9)  is  frequently  used  in  the  literature 
to  obtain  an  approximate  solution  for  an  adiabatic  flow  about 
a  cone  in  supersonic  or  hypersonic  main  stream. 

Joining  of  the  Outer  and  Inner  Solutions 

In  view  of  the  restrictions  of  Eqs.  (7)  and  (8) 
it  is  appropriate  to  Join  the  inner  and  outer  solutions  at 
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Substituting  Eq.  (12)  in  Eq.  (5),  there  is  obtained 


(13) 


In  general  Eq.  (13)  must  be  solved  for  cp,p  by  trial  and  error. 

If  terms  containing  powers  of  (cp  -  cp*)  greater  than  unity 

S  1 

are  negligible  an  analytic  solution  is  obtained. 

Once  (p.p  is  known,  the  coefficients  A  and  B  are  found 
by  equating  Eqs.  (4)  and  (10)  and  Eqs.  (5)  and  (11)  at  cp  -  (p,p 
and  solving  the  resulting  equations.  This  procedure  yields 


B  >■  sin 


(p,j,  v*(cp,j,;q)g,u* 


)  cos  cp,j,  -  u  ((p,j,;(Pg,Ug)  sin  (p,j, 


(1/2)  cos  cp^  log 


1  +  cos 

1  -  cos  (p,j,  y  ^ 


cos  cp- 


(14) 


(15) 


where  u  ((p„;  ^  (‘P'pJ  obtained  from 

X  S  S  X  S  5 

Eqs.  (4)  and  (5). 

From  the  condition  Vg  =■  0,  and  Eq.  (11),  the  relation 
for  determining  the  cone  angle  is  found  to  be 


1  +  cos  (J) 
1  -  cos  cp 


cot  (Pg 
sin  (Pg 


(16) 
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Equation  (16)  must  be  solved  by  trial  and  error  for  the  value 
of  cpg  corresponding  to  a  given  value  of  (-A/B) .  The  plot 
in  Fig.  4  is  convenient  for  ascertaining  a  first  approxi¬ 
mation  to  (Pg,  from  which  a  more  accurate  value  may  be 
obtained,  if  so  desired,  by  trial  and  error. 

COMPARISON  OF  AI9R0XIMATE  AND  "EXACT" 

UNIT  CONICAL  FLOW  FIELDS 

It  is  desirable,  of  course,  to  verify  the  accuracy 
of  the  approximate  unit  field  solution  by  comparing  it  with 
an  exact  solution.  Apparently  the  only  exact  calculations 

performed  for  the  complete  unit  field  are  those  reported  in 

1  2 
Ref.  25  ,  Unfortunately  Ref.  25  records  (p  without  giving 

s 

*  ■!< 

the  corresponding  values  of  u  or  v  required  for  the  present 

s  s 

3 

purpose  .  Examination  of  Table  VI  in  Ref.  25,  however, 
reveals  that  (tpg  -  qi^)  is  generally  less  than  0.2  radians 
for  a  wide  range  of  parameters.  This  suggests  that  u*  can 


The  investigators  in  Ref.  25  performed  these  calculations 
principally  as  a  matter  of  mathematical  Interest  since  they 
considered  the  singular  discontinuity  to  be  without  physical 
significance . 

2 

See  Table  VI  and  Diagram  No.  8  in  Ref.  25. 

3 

An  inquiry  to  MIT  regarding  the  possibility  of  obtaining  the 

values  of  u  or  v  yielded  negative  results, 
s  s 
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be  developed  in  a  series  expansion  In  powers  of  (q>  -  9^) 

from  which  u  can  be  determined  from  the  known  value  of  (|)_ 
s  s 

given  In  Ref.  25.  The  appropriate  expressions  for  this 
purpose  are  derived  In  Appendix  D.  This  procedure  should 

i(i 

give  a  reasonably  good  value  for  u  .  However,  derivatives 

B 

of  u  will  be  less  accurate  near  as  a  consequence  of  the 

singular  behavior  at  9  which  Is  not  accounted  for  In  the 

s 

series  expansion.  For  this  reason  v^  Is  calculated  from 

Eq.  (3)^  rather  than  from  the  derivative  of  the  series  for 
♦ 

u  . 

Exact  unit-field  solutions  also  can  be  obtained  by 
the  more  laborious  process  of  extending  the  flow  fields 
of  Ref.  25  forward  of  the  shock  position  by  means  of  a 
numerical  Integration  procedure.  Using  the  Runge-Kutta 
Integration  technique,  one  such  calculation  has  been  per¬ 
formed  as  a  check  on  the  accuracy  of  the  series  expansion 
method.  In  order  to  provide  a  relatively  severe  test,  an 
example  was  selected  for  which  (q>  -  (Po)  Is  moderately  large 

S  0 

The  parameters  of  the  example  are  (Pg  -  40  deg.,  u^  -  0.40, 

cPg  -  59.431  deg.  and  tp^  -  68.555  deg.  (corresponding  to  the 

flow  field  on  p.  429  of  Ref.  25).  Values  of  q>  «  68.556  deg 

s 

and  Ug  ■■  0.30969  were  obtained  from  the  Runge-Kutta  Inte¬ 
gration.  The  series  expansion  method  gave  a  value  of 
% 

u„  -  0.31050  for  the  retention  of  terms  In  the  series  to 
s 

1  *  ♦ 

Recalling  that  a  -  v  . 

SB  ■>. 
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2 

order  (cp  -  (Pg)  ,  and  a  value  of  -  0.31017  for  retention 

3 

of  terms  to  order  (q>  -  <Pg)  .  It  can  be  seen  from  these 

comparisons  that  the  series  expansion  method  predicts 

Ug  rather  well.  It  can  be  anticipated  that  in  general  the 

agreement  will  be  better  than  in  the  example  for  smaller 

values  of  ((p^  -  <Pg)  and  poorer  for  larger  values. 

"Exact"  values  of  u*  used  in  subsequent  comparisons 

2 

are  based  on  the  retention  of  terms  to  the  order  of  (qp  -  qpg) 
in  the  series  expansion. 

Approximate  and  exact  calculations  of  the  flow 
field  for  the  foregoing  example  are  shown  in  Fig.  5.  In 
the  approximate  calculation  the  outer  solution  was  calculated 
retaining  terms  to  the  order  of  (<p  -  cp)  in  u  .  It  is 

5 

apparent  that  the  approximate  calculation  agrees  very  well 
with  the  exact  one. 

The  results  of  additional  computations  comparing 
exact  and  approximate  values  of  the  cone  semi-apex  angle  and 
the  normalized  cone-surface  velocity  for  given  values  of  <p 

s 

* 

and  Ug  are  recorded  in  Table  II.  Also  given  in  the  table 
are,  the  Joining  angle  cp,p  for  the  inner  and  outer  solutions, 
and  the  free-stream  Mach  number  M  for  which  the  unit  field 

J  1 

corresponds  to  a  Chapman- Jouguet  detonation  .  Note  that  the 


The  free-stream  Mach  number,  M.  ,  is  found  by  substituting 

°J 

the  C-J  condition  in  Eq.  (29)  and  solving  the  resulting 
quadratic  equation  in  M  2. 
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cone-angle  error  (<Pg  -  <I)g)/<l>g  generally  Increases  with  in¬ 
creasing  values  of  ((p„  -  <Po)/cp„,  and  that  the  error  is 
generally  larger  for  the  smallest  cone  angle.  Note  also, 
that  for  a  given  cone  angle,  the  cone-angle  error  decreases 
with  Increasing  free-stream  Chapman- Joug^et  Mach  number. 

The  error  in  the  cone-surface  velocity  is  seen  to  be  very 
small . 

The  cone-angle  error  is  within  acceptable  limits 
for  the  30-  and  40-degree  cones  and  the  10-degree  cone  at 
C-J  Mach  numbers  greater  than  3.0,  but  is,  perhaps,  of 
marginal  acceptability  for  the  10-degree  cone  at  C-J  Mach 
numbers  of  less  than  3.0. 

The  larger  cone  angles  are  of  greater  interest 
for  detonation  experiments  since  their  stronger  shocks  are 
more  conducive  to  producing  a  detonation. 

Approximate  and  exact  calculations  of  the  flow 
field  for  Example  VI  in  Table  II  are  shown  in  Fig.  6. 

THE  DETONATION  FLOW  FIELD 


The  flow  field  under  consideration  is  shown  in 

Fig.  Ic,  with  the  subscript  "0”  denoting  free-stream 

conditions  and  the  subscript  "2”  denoting  conditions  in 

the  field  between  the  exothermic  discontinuity  and  the  cone. 

As  previously  noted,  it  is  assumed  that  M^,  c  , 

'’o 

c  ,  c  ,  c  ,  and  tp ,  are  known . 

^o  ^2  ^2  *' 
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Consider  first  the  determination  of  the  tangential 
velocity  component  downstream  of  the  Chapman- Jouguet  wave. 
From  the  condition  of  constancy  of  the  tangential  velocity 
component  across  the  wave  there  is  obtained  the  relation 


M  (cos  cp_)(a  Vc-) 

O  J  O  ^  j 


where  (a  /c»)  may  be  shown  to  be  given  by 
°  ^  J 


Y  ® 

*  o  o 

Yo 


with  (t  /to)  obtained  from  Eq.  (A16)  as 
o  ^  j 


(17) 


(18) 


(19) 


(20) 
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(21) 


Hiti  al|i 

With  Uj  and  cPj  known,  «Pg  and  may  be  determined 
by  means  of  the  unit  field  solution. 

By  manipulating  Eq.  (3),  which  is  valid  throughout 
field  2,  the  surface  Mach  number  is  easily  shown  to  be  given 
by 


(22) 


The  corresponding  surface  pressure  is 


fB 


1 

(p/p). 


J  ( 


(23) 


where  (p/P)^  is  given  by  Eq.  (A19)  with  M^  replacing  M,  and 
(P/P  )  is  obtained  by  specializing  Eq.  (A17)  to  the  C-J 


condition  and  replacing  "P  and  p  by  P^  and  p^  respectively. 

From  the  aforementioned  unit-field  solution,  values 
of  Uq  and  v^  corresponding  to  an  arbitrarily  assumed  value  of 
cpQ  are  readily  found.  The  corresponding  values  of  and 
(Q/c  T  )  are  then  determined  as  follows.  From  the  condition 

of  constancy  of  the  tangential  velocity  component  across  the 
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detonation  wave,  and  from  relations  for  a  non-Chapman-Jouguet 
detonation  analogous  to  those  of  Eqs.  (18),  (19)  and  (90),  it 
may  be  shown  that 


To  Mq  sin^  <p^ 


U  \ 


2  2 
CSC  (Pjj  -  cot  (Pjj 


/-  \2 
*0 


V  Po 


+  1 


-  0 


(24) 


where,  from  Eq.  (AS) 


-  'l  I'd  “o 


M?  sln^  cp,.  +  f„  -  1 


(25) 


with 


(26) 


(27) 


(28) 
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Solving  Eqs.  (24)  and  (25)  for  1^,  there  is  obtained 


**o  “  ^2 


cot  (p, 


^O  1 


(u/Cg) 


5  -  'Pd  “ 


-1 


The  corresponding  heat  release  is  found  by  substituting 
Eq.  (29)  in  Eq.  (25).  Obviously  cpj^  must  be  restricted  to 

'Pj—  *Pd  ^  “Ps* 

The  surface  Mach  number  as  given  by  Eq.  (22)  is 
unchanged.  The  surface  pressure  is 


Po 


1  /  p  \ 


(29) 


(30) 


Where  (p/P)^  is  given  by  Eq.  (A19)  with  M^  replacing  M,  and 

(P/P  )  is  given  by  Eq.  (A17)  applied  to  a  detonation,  with 
°  D 

P  and  p  replacing  P  and  p. 

The  error  in  the  approximate  detonation  field 
calculation  is,  of  course,  of  the  same  order  as  the  unit- 
field  approximation  upon  which  it  is  based. 


THE  DETONATION-SHOCK  FLOW  FIELD 


It  is  not  the  purpose  here  to  treat  this  flow  field 
(Fig.  Id)  in  any  detail.  The  initial  portion  of  the  expanding 
flow  downstream  of  the  Chapman- Jouguet  wave  may  be  obtained  by 
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selecting  the  plus  sign  In  the  expression  for  In  the  outer 
solution  of  the  unit  field. 

THE  SHOCK-DEFLAGRATION  FLOW  FIELD 

The  treatment  of  this  field  requires  approximations 
for  the  adiabatic  portion  of  the  field,  region  1  in  Fig.  lb, 
and  the  dlabatlc  portion,  region  2  In  Fig.  lb.  There  are  a 
wide  variety  of  approximations  available  In  the  literature 
(see,  for  example,  Refs.  26  through  31)  for  adiabatic  cone 
flow.  It  Is  worth  noting  that  for  most  of  these  the  accuracy 
of  the  approximation  increases  with  Increasing  free-stream 
Mach  number  (see  Ref.  30).  The  unit-field  solution  presented 
herein  also  may  be  used  for  calculation  of  the  adiabatic 
field. 

If  the  shock-deflagration  flow  is  confined  to 
situations  In  which  the  normal -component  Mach  number  ahead 
of  the  flame  is  of  the  order  of  that  occurring  in  turbulent 
and  laminar  flames  (say  <  0.10)  the  dlabatlc  portion  of  the 
field  in  Fig.  lb  may  be  approximated  by  the  Legendre-type 
solution  used  In  the  Inner  portion  of  the  unit  field. 

For  known  velocity  components  just  downstream  of 
the  deflagration,  a  Legendre-type  solution  for  the  region  2 
In  Fig.  lb  yields 

u*  -  A  cos  cp  +  (B/2)  cos  <p  log  t  ~y)  "  ® 


(31) 
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-4> 

V 


A  sin  cp  +  (S/2)  sin  9  log 


1  +  cos  f 
1  -  cos  q) 


6  cot 


9 


(32) 


where 


S  -  sin  <Pp  (Vp  cos  cpp  -  Up  sin  cpp) 


m* 


1  - 


uj  -  B 


/  1  +  cos  <p-  \ 

(1/2)  cos  Oj,  log  I  1  -  ebs-o- j-  1 


cos  CPi 


(33) 


(34) 


The  cone  angle  corresponding  to  a  specified  deflagra¬ 
tion  (hence,  known  values  of  J.  and  §)  may  be  found  by  means 
of  Eq.  (16)  with  A  and  B  replacing  A  and  B. 

The  surface  Mach  number  is  given  by  Eq.  (22). 

The  corresponding  surface  pressure  is 


wft; 


(35) 


where  (P/P  )  is  obtained  from  oblique  shock  relations  (see 
°  S 

e.g.  Ref.  32),  (p/P^)  Is  given  by  Eq.  (A19)  with  M^  replacing 
M,  and  (P/P)p  is  given  by  Eq.  (A17)  applied  to  a  deflagration. 

SOME  PARAMETRIC  CURVES 

A  limited  number  of  parametric  curves  of  the 
properties  of  some  supersonic  flows  with  oblique  discontinuities 
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are  presented  herein  In  Figs.  7  through  12.  These  curves  are 
not  Intended  to  be  comprehensive;  such  a  goal  Is  beyond  the 
scope  of  this  work.  Nor  are  they  Intended  to  be  the  direct 
result  of  applying  the  foregoing  approximations,  since,  where 
possible,  existing  tabulated  fields  are  used  In  part.  The 
purpose  Is  simply  one  of  Illustrating  some  trends  and  magnitudes 
for  parametric  values  within  an  approximate  range  of  possible 
Interest . 

Figures  7  and  8  present  curves  which  are  plots  of 
Eqs.  (AS),  (A8) ,  (A12)  and  (A14),  or  combinations  thereof. 

These  curves,  along  with  the  remaining  equations  of  the 
collection  (A9)  through  (A19) ,  permit  the  determination  of 
the  flow  properties  across  oblique  deflagrations  for 
normal -component  Mach  numbers  ranging  In  value  from  0.0 
to  0.10,  at  a  constant  specific  heat  ratio  of  1.4 

Some  flow  properties  for  the  shock-deflagration 
regime  Illustrated  In  Fig.  lb  are  given  In  Figs.  9  and  10 
for  cone  flows  and  In  Figs.  11  and  12  for  wedge  flows  for 
a  constant  specific  heat  ratio  of  1.4,  free-stream  Mach 
numbers  of  2.0,  4.0,  and  6.0,  flame  normal-component  Mach 
numbers  of  0.040  and  0.80,  and  identical  adiabatic  shock 
angles  for  the  cone  and  wedge  cases. 

Curves  of  the  body  semi-apex  angle  as  a  function 
of  the  normal-component  velocity  ratio  across  the  flame,  for 
constant  values  of  specific  heat  ratio,  free-stream  Mach 
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number,  shock  angle,  and  flame  normal-component  Mach  number 
are  shown  in  Figs.  9  and  11  respectively  for  cones  and 
wedges.  The  flame  Inclination  angles  are  also  Indicated 
on  the  figures.  Shown  in  Figs.  10  and  11  are  the  associated 
cone-  and  wedge-surface  Mach  numbers  from  which  the  surface 
pressure  may  be  calculated  by  means  of  Eq.  (35).  The  cor¬ 
responding  heat  release  parameter  Is  ascertainable  from 
Fig.  7. 

The  adiabatic  portion  of  the  cone  shock-deflagration 
flow  field  required  in  the  construction  of  the  curves  In  Figs. 
9  through  12  was  obtained  from  the  tables  of  Ref.  33.  The 
flame  Inclination  corresponding  to  a  preassigned  stream  Mach 
number,  shock-wave  angle,  and  flame  normal-component  Mach 
number  was  found  by  Interpolation  In  the  table  of  Ref.  33. 

The  cone  angle  corresponding  to  a  given  normal-component 
velocity  ratio  across  the  flame  (Fig.  9)  was  then  determined 
by  means  of  a  Legendre-type  solution. 

The  construction  of  the  wedge-flow  curves  in 
Figs.  11  and  12  involved  a  rather  straight  forward  applica¬ 
tion  of  the  oblique  shock  relations  (Ref.  32)  and  the 
appropriate  equations  in  Appendix  A. 

An  important  point  worth  noting  with  regard  to 
Figs.  10  and  12  is  the  fairly  large  portion  of  the  parametric 
spectrum  for  which  the  body-surface  Mach  number  is  subsonic. 
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This  indicates  that  the  conical  flow  downstream  of  the  de¬ 
flagration  is  partially  subsonic  for  a  fairly  wide  range  of 
free-stream  Mach  numbers  and  cone  angles.  As  previously 
noted,  such  a  flow  cannot  be  attained  in  the  laboratory 
since  on  a  finite-length  model  the  conical  flow  is  destroyed 
by  the  upstream  propagation  of  downstream  disturbances  in 
the  subsonic  field.  In  experiments  aimed  at  attaining  a 
shock-deflagration  type  of  flow,  it  is  therefore  necessary 
to  constrain  the  parameters  to  values  for  which  the  flow  is 
completely  supersonic  downstream  of  the  deflagration.  In 
this  regard  it  is  observed  from  Figs.  10  and  12  that  the 
parametric  spectrum  for  which  the  flow  downstream  of  the 
deflagration  is  completely  supersonic  Increases  in  range 
with  increasing  free-stream  Mach  number. 

For  cone  flow,  curves  of  the  Chapman- Jouguet  wave 
angle  as  a  function  of  free-stream  Mach  number  for  various 
cone  angles  and  a  constant  specific-heat  ratio  of  1.405 
are  shown  in  Fig.  13  along  with  the  corresponding  adiabatic 
shock-wave  angles.  The  detonation  curves  were  obtained  by 
extending  the  flow  fields  given  in  Ref.  25  forward  to  the 
known  singular  discontinuity  surface  of  Ref.  25  by  means 
of  the  series  expansion  method  described  in  the  section 
"Comparison  of  Approximate  and  Exact  Unit  Conical  Flow 
Fields."  The  previously  described  similarities  in  the 
characteristics  of  the  detonation  and  adiabatic  shock  waves 
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are  readily  apparent  in  the  figure.  The  condition  for 
which  sonic  velocity  occurs  at  the  cone  surface  Is  Indi¬ 
cated  by  the  circular  and  triangular  symbols  respectively 
for  detonation  and  adiabatic  shock  waves.  The  corresponding 

cone  angles  are  opp  and  cp.  .  As  noted  on  the  figure,  the 
“j  o 

surface  Mach  number  Is  subsonic  for  detonation  wave  angles 
greater  than  the  wave  angle  for  a  sonic  surface  velocity, 
and  supersonic  for  lesser  wave  angles. 

Since  a  sonic  surface  velocity  Is  the  limiting 
condition  for  the  attainment  of  conical  flow  on  a  finite 
length  model  In  an  experimental  facility,  an  "a  priori" 
knowledge  of  qp-  Is  Important  to  the  selection  of  a  test 

model  for  which  the  semi-apex  angle  should  be  as  large 

as  possible  to  promote  detonation  without  destroying  the 

conical  flow.  The  adiabatic  cone  angle  op^  Is  not  satls- 

“o 

factory  for  this  purpose  since  It  Is  generally  larger  than 

cPq  .  Note  on  Fig.  13,  for  example,  that  (p^  30  degrees 

®J 

for  a  free-stream  Mach  number  of  2.2,  whereas  the  correspond- 
Ing  adlabatlc-f low  cone  angle  cp-  at  the  same  Mach  number 

Is  found  to  be  39.1  degrees  (Ref.  32). 

It  is  well  known  for  adiabatic  cone  flow  that 
there  exists  an  absolute  maximum  cone  angle  above  which  a 
conical  flow  region  cannot  exist  regardless  of  the  free-stream 
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Mach  number.  For  a  specific  heat  ratio  of  1.405  the  value 
of  this  maximum  angle  is  known  to  be  57.6  degrees  (Ref. 

25) .  The  analogous  absolute  maximum  cone  angle  for  a 
C-J  detonation  flow  was  found  herein  to  be  35.8  (-  0.2) 
degrees . 


CONCLUDING  REMARKS 

Approximate  relations  have  been  developed  for 
calculating  cone  flow  fields  of  the  shock-deflagration  and 
detonation  types  with  attached  adiabatic  and  diabatic  dis¬ 
continuities. 

The  appropriateness  of  the  Lengendre-type  approxima¬ 
tion  used  for  the  shock-deflagration  flow  has  been  demonstrated 
various  places  In  the  literature  in  applications  to  adiabatic 
cone  flow.  The  appropriateness  of  the  present  approximate 
method  for  detonatlve  flow  fields  has  been  verified  herein 
where  it  is  shown  that  the  error  decreases  with  increasing 
cone  angle  and  Increasing  free-stream  Chapman- Jouguet  Mach 
number. 

It  is  believed  that  the  analytical  relations  of 
this  paper  will  be  useful  in  the  design  and  interpretation 
of  cone-flow  supersonic  combustion  experiments  since  there 
are  no  flow  field  tables  available  for  the  specific-heat 
ratios  and  conditions  likely  to  be  encountered  in  experiment, 
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and  hand  numerical  Integration  of  the  governing  differential 
equation  Involves  a  greater  effort  than  does  the  present 
method.  The  availability  to  the  researcher  of  a  digital 
computer,  of  course,  circumvents  the  use  of  the  present 
method.  Such  facilities  are  not  always  available,  however, 
and  in  any  event  the  present  method  serves  as  a  useful 
complement  to  machine  calculation  as  well  as  providing  use> 
ful  relations  for  the  machine  program  in  the  vicinity  of 
the  singular  discontinuity. 
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APPENDIX  A 

Single  Oblique  Exothermic  Discontinuity 

Consider  the  deflagration  shown  in  Fig.  2  to  be 
any  general  exothermic  dlscontipuity  (detonation  or  de¬ 
flagration)  separating  two  conical  flow  fields.  In  view 
of  the  assumption  of  a  perfect  gas  in  the  separate  regions, 

the  following  thermodynamic  relations  and  definitions  apply 

2  . 

in  a  given  region:  p  -  pRT,  R  -  Cp  -  c^,  a  -  y'RT  -  yp/p, 
and  Y  -  Cp/c^ . 

With  reference  to  Fig.  2,  the  equations  of  con¬ 
servation  of  mass,  normal  momentum,  tangential  momentum,  and 
energy  across  the  dlscont-lnuity  (denoted  by  the  subscript  d) 
are  respectively. 


Pd  ^d 


- 


Pd  ^d 


(Al) 


Pd  +  Pd  '^d 


m  m  m2 

Pd  +  Pd  ^d 


(A2) 


"d  -  '‘d 


(A3) 


1  _2  1  .2 
o  V .  +  c^  T  +  Q  -  V .  +  C^  T . 
2d  pj  d  2d  pg  d 


(A4) 
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The  foregoing  derived  forms  of  the  momentum  and  energy 
equations  are  obtained  from  the  basic  forms  through  the  use 
of  the  conservation  of  mass  relation  In  the  basic  momentum 
equations  and  the  conservation  of  mass  and  tangential  momentum 
In  the  basic  energy  equation. 

The  equation  of  state,  the  relations  m^ 

and  "  T  R/(r  -  1)  and  Bqs.  (Al),  (A2)  and  (A4)  may  be 
combined  to  give  t 


'd 


(A5) 


is 


Equation  (A5)  may  be  solved  for  (v/v)^.  The  result 


V 

V 


^7^  ^2  “d  y  ± 

irz  +  1)  ®d 


+  F 


(A6a) 


where 


F  - 


(^2  -  rjXrg  +  1) 

L  ''i  '>'1  ■ _ 

_2 

“d 


_2 
-  2  ®d 


(A6b) 
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and 


“d  "  "“d  ^^d  "  "^d^ 


(A7) 


Discussion  of  the  consequences  of  the  selection  of 
a  specific  sign  preceding  the  radical  In  Ec.  (A6a)  Is  relegated 
to  Appendix  B. 

Equation  (A3)  and  the  velocity  vector  diagram,  Fig. 

2,  yield  the  following  result  for  the  change  In  flow  direction 
across  the  discontinuity 

tan  (cp^  -  e^)/tan  (cp^  - 


Equations  (A6)  and  (AS)  are  the  governing  relations  for  an 
oblique  dlabatlc  discontinuity.  For  known  values  of  m^, 

Q/Cp  T^»  and  t**®  foregoing  equations  may  be  solved 


for  (v/v)^  and  0^.  In  some  problems,  different  combinations 
of  these  parameters  may  be  specified. 

The  density,  static-pressure,  static- temperature, 
sonic  speed,  and  normal  Mach-number  ratios  across  the  dis¬ 
continuity  may  be  expressed  in  terms  of  (v/v)  ..  Relations 

d 

for  the  preceding  quantities,  derivable  from  Eqs.  (Al)  and 
(A2)  and  the  basic  thermodynamic  relations  and  definitions  are 
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Co/p)j  - 

(A9) 

(p/p)d  - 1  Il"^^d "  3 

(AlO) 

(T/T)^  -  (Rj/R2)(v/v)^(p/p)^ 

(All) 

F 

(a/a)^  -  (y2/ri)(v/v)^(p/p)  J 

(A12) 

1 —  — 1 

(m^)^  -  (^/^) ^/ (p/p) ^ 

(A13) 

where  (p/p)jj  in  Eqs.  (All),  (A12) ,  (A13),  is  given  by  Eq.  (AlO) . 

The  resultant-velocity  ratio,  the  resultant  Mach- 
number  ratio,  the  stagnation- temperature  ratio,  and  the 
stagnation-pressure  ratio  depend  upon  vector  quantities. 
Utilization  of  Fig.  2,  Eqs.  (A3)  and  (A4)  and  basic  defini¬ 
tions  yield 


(V/V)^  “  "  ^)jj/cos  (cp^  -  e)^  (A14) 

(M/M)^  -  (V/V)^(¥/a)^  (A15) 
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(5/P)^  -  (p/p)^  (p/P)d/(p/P)ci 


where 

(T/t) 

(p/P) 


(A17) 

(A18) 

(A19) 


yield  the  required  quantities  In  Eqs.  (A16)  and  (A17)  by 

appropriate  substitution  of  and  for  M. 

A  wide  variety  of  alternate  forms  of  the  preceding 

relations  may  be  obtained  by  suitable  manipulations. 

For  specified  values  of  y, ,  Yn>  and  Q/c  T., 

X  z  Pi  n 


the  flow  quantity  ratios  given  by  Eqs.  (A6)  and  Eqs.  (A9) 
through  (A13)  are  functions  of  the  normal  component  of  the 
upstream  Mach  number  in  a  manner  analogous  to  the  same  well- 
known  property  for  an  oblique  shock  wave  (adiabatic  dis¬ 
continuity)  .  The  stagnation  pressure  ratio,  however,  does 
not  follow  this  analogous  behavior  in  the  diabatic  case. 

Equations  (A6),  (AlO) ,  and  (A13)  are  analyzed  in 
Appendix  B  for  the  purposes  of  delineating  the  boundaries  of 
the  various  normal-component  flow  regimes  and  classifying 
them  according  to  Jouguet's  rule. 
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APPENDIX  B 


Delineation  of  the  Normal-Component  Flow  Regimes  for 
an  Oblique  Exothermic  Discontinuity 


For  the  discussion  in  this  appendix  it  is  convenient 
to  rewrite  Eq.  (A6),  introducing  a  special  notation,  as  follows 


where 


>2 


_2\ 


yT  +  rg  «d]±  ^ 


+  1) 


(Bla) 


ffl  .  .  2 

(^2  ■  ^1^^^2 

1)  -j 

F  - 


d 


-  2 

a 


>2  - 


Q 


“•’A 


(Bib) 


and  the  sign  of  the  quantity  (1)  in  the  subscript  of  (v/v)^ 
is  taken  to  agree  with  the  sign  selected  ahead  of  the  radical. 
Quantities  which  are  functions  of  (v/v)^  also  carry  the  same 
subscript  notation. 

From  Eq.  (Bl)  it  is  seen  that  (v/v)^  is  double 

valued  for  given  positive  values  of  F  and  m^,  that  there  are 

two  possible  real  values  of  corresponding  to  given  positive 

values  of  F  and  Q/c  T^,  and  there  are  no  real  solutions  for 

Fr*> 
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F<  0.  It  Is  apparent  therefore  that  the  factors  crucial 
to  the  delineation  of  the  various  normal-component  flow 
regimes  by  means  of  Eq.  (Bl)  are  the  magnitude  of  the 
function  F  and  the  choice  of  sign  preceding  the  radical. 

For  F  -  0,  Eq.  (Bla)  becomes 


/v\ 

vIj  ^^2  +  “d 


(B2) 


Substitution  of  Eq .  (B2)  In  Eqs.  (AlO)  and  (A13)  results  In 

m^  ■■  1.0;  the  Chapman- Jouguet  condition. 

From  the  expression  for  the  derivative  of  (v/v)^ 

with  respect  to  It  Is  easy  to  show  that  the  C-J  condition 

(F  -  0)  Is  a  singular  point  of  Eq.  (Bl)  for  a  constant  Q. 

It  Is  also  easy  to  show,  by  means  of  Eq.  (AS),  that  the  C-J 

condition  corresponds  to  the  maximum  heat  release  for  a  given 

value  of  (v/v)  ,. 

d 

The  heat  release  corresponding  to  the  Chapman- Jouguet 
condition  is  given  by 


where 


C 


r 


^2 

ri  *  (>-1 


(B3) 


(B4) 
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Or,  alternatively,  the  normal  Mach-number  component  yielding 
the  Chapman-Jouguet  condition  for  a  given  heat  release  is 


-2 


"■r[±] 


-  (Cy  c,)  ± 


(Cy  +  C,)‘  -  (rj/r 


.'■1 


(B5) 


where 


(B6) 


and  the  sign  of  the  quantity  L— Jin  the  subscript  of  mj  is 
taken  to  agree  with  the  sign  used  ahead  of  the  radical. 
Quantities  which  are  functions  of  nij  also  carry  the  same 
subscript  notation.  Note  that  the  parenthetical  subscript 
notation  (i)  pertains  to  Eq.  (Bl)  and  quantities  derived 
therefrom,  whereas  the  bracketed  subscript  notation  per¬ 
tains  to  Eq.  (B5)  and  correspondingly  derived  quantities. 

The  effect  of  sign  selection  in  Eq.  (Bl)  upon  the  normal 
component  of  the  downstream  Mach  number  is  now  examined. 
Since  (p/p)^  is  a  positive  quantity*,  and  > 

(v/v)d(_),  it  follows  from  Eq.  (AlO)  that  (p/p)jj^_^j 
(p/p)^j(_)-  Eq.  (A13)  then  yields  m^  (+)  ^  “d(-)’ 

C-J  condition  is  the  boundary  between  the  plus  and  minus 
solutions  of  Eq.  (Bla)  there  is  obtained 

“d(-)  ^  i"d(+) 


(B7) 


wm  «  /  1 

The  ratio  (p/p)^  is  positive  for  (v/v)^  <  1  +  - 

‘'l  ”d 
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That  is,  selection  of  the  minus  (-)  sign  In  Eq ,  (Bla)  gives 
a  subsonic  normal  component  for  the  downstream  Mach  number, 
whereas  selection  of  the  plus  (•»•)  sign  gives  a  supersonic 
normal  component  for  the  downstream  Mach  number. 

Two  pairs  of  curves  given  by  Eq.  (Bl)  for  constant 
values  of  Q,  say  and  >  Q^)  are  Illustrated  In 

Fig.  3  for  >  >'2 *  locus  of  the  C-J  points  given  by 

Eq.  (B2)  Is  shown  by  the  curve  ABC.  It  Is  obvious  (as  noted 
on  the  figure)  that  regions  above  the  C-J  curve  correspond 
to  the  selection  of  the  plus  (-f)  sign  In  Eq.  (Bla)  and  the 
regions  below  the  curve  correspond  to  the  selection  of  the 
minus  (-)  sign. 

As  previously  noted  there  are  no  real  solutions 
to  Eq.  (Bla)  for  F  <  0.  Applying  the  condition  P  <  0  to 
Eq .  (Bib),  the  normal-component  Mach  number  range  for  which 
there  are  no  solutions  to  Eq.  (Bla)  for  specified  values 


of  ^2’  ^d 


Is  found  to  be 


-2  .  -2  ^  -2 

"d  <  ”j[+] 

where  given  by  Eq.  (B3).  The  “jJ+'j 

dltlons  for  Q  -  and  occur  at  the  points  J'  and  J ' 
respectively  In  Fig.  3.  The  corresponding  nij conditions 
occur  at  the  points  J  and  j  In  the  figure. 
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From  Eqs.  (B7)  and  (B8)  It  Is  apparent  that  there 

exists  a  heat-release  quantity  at  and  below  which  real 

solutions  to  Eq.  (Bl)  are  obtained  for  all  real  values  of 

This  occurs  when  m,r-i  "  fo*"  which  the  term  under 

d  j[-]  j[+] 

the  radical  in  Eq.  (BS)  becomes  zero.  The  corresponding  values 
(Q/c  m-^®\  and  (v/v)  are  easily  shown  to  be 

Pi  u  J  J 


(0) 


m 


J 


1 

(?'2/>'l>^ 


(v/v)j^®^  «  1.0 

The  velocity  ratio  curves  are  given  by 


1.0 


I ^yo)  ^  2  rg/ri  (rg  -  D 


> 


-  (0) 
"»d  ^  "J 


(B9) 


(BIO) 

(BID 


(B12) 
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2 

2  +  (^2  -  1) 


1.0 


I  -  _(0) 

'  “d  ^  “j 


(B13) 


Curves  for  Q  -  are  shown  in  Fig.  3.  Point  B 

in  the  figure  corresponds  to  . 

If  the  parameters 

unity,  the  following  asymptotic  forms  for  Eq.  (Bl)  are  readily 
obtained . 


md->0 


r2  + 


^  ^1 


m 


(Cy  +  C^) 


1 

j 

^d(-) 

mq-^O 

*  7 

^  (Yz/Y^)^ 

V 

V 


d(+) 


-  1.0 


"d'^* 


/ 


(-) 


— >  “> 


^2  -  ^ 


(B14) 


(B15) 


(B16) 


(B17) 
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Finally,  the  limiting  boundaries  for  the  regimes  of 

exothermic  discontinuities  are  delineated  by  the  adiabatic 

condition  Q  •-  0.  For  the  adiabatic  condition,  note  from 

Eq.  (Bl)  that  F  is  always  positive  and  greater  than  zero, 

providing  Therefore,  for  Q  -  0  and  real 

solutions  to  Eq.  (Bl)  are  obtained  for  all  real  values  of  m.. 

d 

In  Fig .  3  the  endothermic  regions  are  Indicated  by 
the  hatched  areas.  The  shaded  areas  delineate  regions  of  weak 
deflagrations  and  detonations,  while  the  unshaded  areas 
represent  strong  deflagrations  and  detonations.  The  bounding 
adiabatic  curves,  Q  0,  are  noted. 

Although  the  material  presented  in  this  appendix 
is  not  new  in  its  entirety,  it  does  cover  the  subject  in  a 
slightly  more  generalized  manner  than  heretofore. 
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APPENDIX  C 
The  Omega  Functions 

Oj  -  I  (1  +  d’ 

CI2  -  -(5)(i  +  2r)/p^ 
ng  -  -(1  -  4r)/Pj 

-  (17  -  lor  +  8r^)/F2 
fig  -  (41  +  32r  -  16r^)/F2 
ng  -  (-58  -  22r  +  8r^)/F2 
Clj  -  (-320  -2,280r  -240r^  +  320r^)/F3 
Og  -  (-5,040  -10,440r  -5,400r^)/F3 
Og  -  (150  +  3,780r  -8,840r^  -960r^)/F3 
n^Q  -  (1,920  -3,240r  +  18,000r^  +  960r^)/F3 
Oil  -  (-1,750  +  l,740r  -9,120r^  -  320r^)/F3 
where 

r  =  (r  -  i)/2  Fj  a  12  (1  +  r) 

F2  =  270  (1  +  r)^  F3  -  72,900  (1  +  r)^ 
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APPENDIX  D 


Denote  the  angular  position  of  an  arbitrary  non- 
singular  conical  ray  by  q)  ,  then 

A 

u*  -  u*  -  V*  (q)  -  q)^)  +  Ag  (q)  -  +  Ag  (q)  -  q)^^)^  + -  (Dl) 

where  (q>  -  q)^)  is  a  small  quantity.  Substitution  of  (Dl)  in 
Eqs.  (1)  through  (3)  yields 


where 

k'  -  cot  q)  . 

A 

The  corresponding  v’*'  component  is  found  by  means 
of  Eq.  (2). 
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Numerical  Values  of  the  Omega  Functions 
for  Several  Specific-Heat  Ratios^ 


Y  -  1.2 

Y  -  1.3 

y  >■  1.4 

Y  “  1.405 

^1 

1 . 5732 

1.6086 

1.6432 

1.6449 

^2 

-0.45455 

-0.47101 

-0.48611 

-0.48683 

-0.045455 

-0.028986 

-0.013889 

-0.013167 

"4 

0.049219 

0.043912 

0.039403 

0.039196 

^5 

0.13480 

0.12726 

0.12027 

0.11993 

^6 

-0.18402 

-0.17117 

-0.15967 

-0.15913 

^7 

-0.0056692 

-0.0060097 

-0.0062160 

-0.0062234 

^8 

-0.063259 

-0.060677 

-0.058299 

-0.058185 

^9 

0.0045413 

0.0046842 

0.0043877 

0.0043640 

^10 

0.018314 

0.016616 

0.015874 

0.015857 

"ll 

-0.017186 

-0.015290 

-0.014046 

-0.013997 

^See 

Appendix  C  for 

omega  function 

relations . 
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Fig.  1  CONICAL  FLOW  REGIMES  FOR  INSTANTANEOUS  CHEMICAL  REACTION 


Fig.  2  ILLUSTRATION  OF  NOTATION 


Fig.  3  SCHEMATIC  ILLUSTRATION  OF  NORMAL-FLOW  REGIMES  FOR  OBLIQUE 
EXOTHERMIC  DISCONTINUITIES 
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Fig.  12  TODGE-SURFACE  MACH  NUMBER  FOR  SHOCK-DEFLAGRATION  FIX>WS 
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Fig.  13  WAVE  ANGLES  OF  CHAPMAN -JOUGUET  DETONATIONS  AND  ADIABATIC 

SHOCKS  FOR  CONE  FLOW 
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